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1/ Introduction 


In the 1987-1990 period NASA has planned several ground and flight 
experiments with the eventual objective of deploying large flexible structures 
in space. A currently active precursor is the SCOLE experiment [1]. Here 
the problem is that of slewing an offset antenna on a long (130 ft.) flexible 
beam-like truss attached to the space shuttle, with rather stringent pointing 
accuracy requirements (±.02 degrees). This paper examines the relevant 
methodology aspects in robust feedback-control design for stability augmen- 
tation of the beam using on-board sensors. We frame it as a stochastic 
control problem - boundary control of a distributed parameter system described 
by partial differential equations. While the framework is mathematical, the 
emphasis is still on an engineering solution. 

The fact that the deployment is in space makes model uncertainty the 
major consideration in control design. Particularly serious in this regard 
is for instance the modelling of inherent damping in the system long known 
to be difficult [2] , and a still unresolved problem even in theory. Hence 
robustness becomes a must feature, even at the expense of optimality. Another 
aspect is the complexity of computation, making any simulation study a costly 
undertaking. 

The overall model involving both slewing and beam stabilization is still 
not well understood. Hence the two problems — of slewing and stabilization — 
are best studied, at least in initial efforts such as reported here, separately. 
We attempt stabilization at the termination of the slewing so that in parti- 
cular the system is essentially linear except for a small nonlinear term 
contributed by the kinematic nonlinearity. It should be noted that at present 
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we do not have a stochastic time-optimal control theory adequate for optimal 
slewing based on sensor data. 

An abstract mathematical formulation is developed in Section 2 as a 
nonlinear wave equation in a Hilbert space. We show that the system is 
controllable and develop a feedback control law that is robust in the sense 
that it does not require quantitative knowledge of system parameters. The 
stochastic control problem that arises in instrumenting this law using 
a PP ro P r iate sensors is treated in Section 3. Using an Engineering first 
approximation which is valid for "small" damping, formulas for optimal choice 
of the control gain are developed. 

2. Abstract Formulation 

We are concerned with the mast stabilization problem only and the model 
we use assumes that the angular velocity of the shuttle-antenna system is 
small enough to be neglected. We model the mast as a thin prismatic beam. 
There is then the question of whether a finite-element model or a continuum 
(involving partial differential equations) model should be used. Here we 
deal only with the latter, the basic governing equations being beam bending 
and torsion equations with controls at the boundaries. 

With reference to Figure 1, the beam of length L is along the Z 
axis, z being zero at the shuttle end. u^(-), u 0 (*) will denote the 

displacements along the Y-Z, x-z planes and u^(*) the angular deflection 
about the Z axis. In addition proof-mass controllers are provided at 
points s 1 and s 2 , on the beam, the locations to be chosen optimally. 
Control moments are applied at both ends as well as control forces at the 
reflector center. The various moments of inertia and masses are specified 
in [1], (2]. 
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FIGURE 1 

SHUTTLE/ANTENNA CONFIGURATION 
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We first develop an abstract mathematical model. We define 
H - L 2 rO,L ) 3 x R 14 o < L < co 

with the usual inner-product thereon denoted [ , ]. We fix the points 

° < s 2 < s 3 < L and define a linear operator A into H with domain D 
m H defined as follows. We use u^(-), u 0 (-), u^-) to denote the 

functions in L_,[0,L] 3 . Thus an element x in H is denoted 

V ,} 

v*> 




X 17 



The domain 

D consists of elements 

x such 

that u', 
<P 

u e' € V°' L] 

and u'" ( • ) 

has - derivatives in 

[0,s 2 ] , 

fS 2 ,S 3 1 

and [s 3 ,L]; 


similarly for u Q ( • ) ; ( • ) such that u^(*) and uj(-) « L 2 [0,L] ; the 

remaining components of x are specified as 


X 4 
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u. (0+) 

V 

X 11 

= 

u;(L-) 

V 

X 5 

= 

u e (0+) 

X 12 
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u e (L_) 

X 6 

= 

u.(l-) 

X 13 
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X 7 

= 

u e ( L-) 

X 14 
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% (s 2 ) 

X 8 
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u! (0+) 
4 1 

x is. 

s 

u e (s 2 ) 

X 9 
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u 0 (0+) 

x l6 

- 

W 

X 10 
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V 0+ > 

x,. 

m 

u»(s.) 


Thus at least for x in D, we may identify the finite-dimensional part 
as the "boundary." The operator A is then defined by 

y * Ax 

where the functional part (in L 2 (0,L] 3 ) is given by 

ei u;-(.) 

9 ♦ 

EI e u e" ( * ) 

-GI^uJC) 

and the boundary part by: 
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X 
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y n 
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y 5 

s 

El gUg" (0+) 
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♦ $ 
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X 

G V* (L " ) 

y 7 

= 

- EI e u e* (L -> 

y i4 

- 

E V u i ,,(s 2 +) - u ; ,,(s 2- n 

GO 

>1 

- 

-El u"(0+j 
9 9 

y i5 

= 

EI e (u e' (s 2 +) - u 9" (s 2‘ )) 

y 9 

- 

" EI e u e (0+> 

y i6 

* , 

E V u r ( v> - u ; ,,(s 3’ )) 

y io 

* 

- G V* (0+) 

y !7 

m 

E1 e (u e” (s 3 +) ‘ u e ,(s 3 - )) 


It may then be verified that D is dense and A is self-adjoint and 
nonnegative definite. Moreover A has a compact resolvent with a complete 
orthonormal set of eigenfunctions (modes). Zero is an eigenvalue. 

The control system dynamics can then be characterized as a nonlinear 
wave-equation: 

Mx(t) + Ax ( t) + K (x ( t) ) + Bu (t) = 0 (2.1) 

where M is a 17 * 17 nonsingular nonnegative definite matrix, and defines 
self-adjoint positive definite linear operator H onto H. The control 

u(*) is in R 12 , and 
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Bu 


X - co1 - (0 ' U 1' u 2 , u 3> u 4 , u 5 , u 6 , u 7 , u 8 , u 9 , u 10< u xl , u 12 I 

We have thus only "boundary" control. The nonlinearity is kinematic: 


where 


K(x) 


0 

0 

0 

0 

0 

0 

0 


n ! • 'A 

n 4 • V. 


0 

0 

0 

0 


n l = 

col 

(x e ,x 9 ,x 

ft 

4 

col 

(X 11' X 12 


^4 are symmetric positive definite (moment) matrices and 
• denotes vector cross-product. 

Two relevant properties of the function K(-) are: 

(i) [ K (x) , x) = 0 

(ii) 1 K ( x ) || < const. ||x|| 2 

We do allow for "state noise" and let 


4 


A 
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A , and the 


N(t) 


FN (t) 


N x (t) 
N 2 (t) 
N 3 (t) 


x ( t) 


where N(t) is white Gaussian with spectral density matrix 
components of x(t) are defined by 


i = 1 f • • • f 7 


x i (t) = 0 

Xg(t) = ^(t) 

x 9 (t > = N 2 (t) 

X 10 (t) " N 3 (t) 

x.(t) =0 i > 10 . 

Note that the "boundary" values are part of the state. 

State-space Form 
With 


Y(t) = 


x(t) 
x (t) 


we go over to the state-space form: 

Y(t) = Ay( t ) + K(Y(t)) + 8u(t) + F(N(t)) 

where 


8u(t) 


-M -1 A 0 


-M _1 Bu(t) 


( 2 . 2 ) 
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and in the notation 


we have 


Y = 



Y c // x // 


ICC Y) 


FN(t) 


0 

-M -1 K(y 2 ) 

0 

-M -1 FN (t) 


well known, we can introduce a new inner product, the "energy" 
inner product 

[Y,Z] £ = [7a y x , ✓a zj + [My 2 , z,,) 

on /?(A) x « . i?(A) is the orthogonal complement of the null space of A. 

We denote the completed space by We shall from now on consider only 

We have: 

A + A* = o 

and of course A has a compact resolvent and we have an orthogonal decompo 
sition of H given by 


Y = 



where P R is a two-dimensional projection for each k, 



A Vk 


( 2 . 3 ) 


spanned by 
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where 




*K • 


2 ^ _ 2 

\ > 0* <*> k •; (2. 


tM* k , ♦ J 


‘ 5 - 


Let S(t) denote the semigroup generated by A. Then we have the 
sentation: 


repre- 


S (t ) Y = 


l S(t) P Y 
1 K 


More explicitly, if 


Then 


P k S(t) P k 


S (t) Y = 


s( t) P R . 


y x (t) 


y 2 (t) 


y 2 (t > = y 2 (t) 


and 


CO 

V l (t) * l fy l' M *k ] *k COS V + [ ty 2' M V * k oT 


sin 


V 


Note that it is required that y satisfy: 


(2.5) 


l ty r M# k ] 0,2 < 


It is easy to establish existence and uniqueness of solution for the 
integral version of (2.2): 

t t 

Y ( t) * S (t) Y (0) + / S(t-o) 8u(o) do + / S(t-o) Fn(o) do 

0 0 

t 

+ / S(t-o) K(y (o ) ) do , (2.6) 

0 

without invoking any nonlinear semigroup theory, by just Picard iteration. 
See [3] . 


377 



We can now state the basic result that yields a robust feedback-control 
law for the deterministic system (seeing F = 0) . 


Theorem 2.1. 


Let P be any 12 * 12 symmetric nonnegative definite nonsingular matrix. 
Then the feedback control 


u(t) = -P 8* Y( t ) 

is such that the "closed-loop" system 

Y(t) = Ay(t) - 8pB*Y(t) + K(Y(t)) 
is globally asymptotically stable. That is to say 


l*<t)fl E + 0 


(2.7) 


( 2 . 8 ) 


as t 00 

Proof. We refer to [4] for a proof. The proof exploits the fact that (A, 8) 
is controllable in an essential way. in particular the 
generated by (A - 8 P 8*) is strongly stable: that is to say: 


semigroup S (t) 

B 


|S B (t)Yfl E ■* 0 


as t + 00 . 


We also obtain that 


l ( p 8 *Vt)Y, 8 *s B (t)Y) dt = | ||y||2 . 

The control law is also optimal for the quadratic cost functional: 

°° OO 

/ I|/F B*Y(t) II 2 dt + / llu(t) II 2 dt 
0 0 


for the linear system 


Y(t) = Ay (t) + 8/p u (t) 


(2.9) 
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3. Stochastic Control 


To instrument the control law 

u (t) - PB*x(t) 

- Pb(t) . (3.1) 

We need to assume co-located (rate) sensors. The sensor output v(t) 
would then be: 

v(t) = b(t) + N Q (t) (3.2) 

where N Q (t) represents the sensor noise, modelled as white Gaussian with 
(12 * 12) spectral density matrix D. In terms of the state-space represen- 
tation (2.2), we can rewrite (3.2) as 

v(t) = CY ( t) + N Q (t) (3.2) 

where 

C = B* 

and C is of course finite-dimensional. If we assume that the separation 
principle applies, a reasonable choice of control law would be 

u(t) a Pb(t) (3.3) 

where, E denoting conditional expectation: 

b(t) = E[b(t) j v ( s) , s £ t) 

and of course 

b(t) = CY ( t) 

a 

where Y(t) is the Kalman state estimates 

Y(t) = E[Y(t) | v (s) , sit] . 
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Even if we were to neglect the nonlinear term K(-), this would require an 
infinite-dimensional Kalman filter, which even if we could instrument it, 
would depend on quantitative knowledge of the system parameters. Hence this 

filter would need to be simplified in considerable measure, in favor of 
robustness . 

The simplest version would be one that did not distort b(t) and thus 
would lead to the control law: 

u(t) = Pv(t) . (3.4) 

We are thus introducing a noise input into the system which may excite 
higher-order modes. Let us therefore study the system response which is 
now given by the stochastic equation: 

Y(t) = (A - 8p8*) y ( t ) - Bptyt) + K(y ( t) ) + Fn ( t) (3.5) 

This can be expressed as an integral equation: 

t 

Y(t) = Y 0 (t) + / S B (t-o) K( Y(o)) do (3.6) 

where 

t t 

Y o (t) = s B (t ) Y < 0 > - / s (t-o)BPN (o)do + / s (t-o) F n(o) do 

0 u 0 B 

(3.7) 

We note that because K( • ) is locally Lipschitzian , we may solve (3.6) by 
Picard iteration: 

t 

Y n +1 = Y o (t) + / S B (t-o)K(Y n (o)) do . (3. 8) 

We omit the details; see [ 3 ]. More important to us is actually (3.7). 

We want to show that the process Y^-) is asymptotically stationary and 
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evaluate its covariance function. Following [5], since S B (*) is strongly 

stable, it is only necessary to show that for Y in Hi 

E 


/ is (0)Bpdp8*S.(0)*Y, Y]_ do < 

o B B E 


and also that 


/ [S (o)FaFV (o)*Y, YJ _ da < 
0 B B E 


(3.10) 


For this purpose we note that Sg(t)* is strongly stable with generator 

A* - 8p8* 


and analogous to (2.9) we have that 


J ll^e*S n (t)*Y|| 2 dt = i |y|| 2 . 


Hence 


Since 


oo 

/ | Sd p8*s (t)*y|| 2 dt J 


< - . 


I f *s b (t) * y| < |8*s B (t)*Y| 


we also obtain (3.10) . For Y, Z in H let 

E 


[R(t,s)Y, Z] = E ( [Y q (t) , Y ] [ Y q ( s ) , Z ] ) . 


Then we have that 


R(t,s) = S (t-s) R(s,s) , 

D 


t £ s 


and bonce it follows that 


limit R (t+L, s+L) « S (t-s) R , t i s (3.11) 

L-*® a " 


* / |/Drt*S (t)*v| 2 dt ♦ / |/AF*S.(t)*V| 2 dt . (3.12) 
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The process Y Q ( * ) is thus asymptotically stationary with covariance operator 

S B (t ' s)R « t * s 

note that is not necessarily nuclear, even though R(t,t) will be if 

R(OfO) is. Indeed taking 


D = dl ; 


P = I 


we obtain that 


/ | /d pB*s b (t) *y| 2 dt = | JyK 2 . 

0 2 E 


From (3.8) we can show that the process Y(t) is asymptotically stationary, 
since y.) will have this property for each n. Since it would appear 
that the nonlinearity is small, we shall now concentrate our attention on 

the linear approximation Y ( • ) . 

o 

The eigenfunctions of (A - BpB*) * ^ 

) are approximately the same as that 

of A and the eigenvalues are 

a k 


°k 1 iw k ' 


« 1 


* where 


Hence 


2 °x - [pb k' b k ' • 


+ + , 


lR,f f, , J (IDPb k' Pb K^I-*,.r-K]) 
k k E k [Pb k' b k ] 


(3.13) 


(3.14) 


which is thus the noise energy in the kth mode. We see that increasing 

the damping but also increases the noise excitation. In prac- 
tice one would want a compromise between increasing damping at selected 
low order modes but keeping the noise excitation ,t higher order modes within 
hound. Clearly further work is needed before any attempt at control design. 
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We may also mention one point of purely theoretical interest. To 
characterize the distributions of the noise response of a nonlinear system 
described by ordinary differential equations one uses the Fokker-Planck- 
Kolmogorov equations which are partial differential equations. In (3.5) 
we have a nonlinear partial differential equation; it would be of interest 
to develop a corresponding tool to study the distributions. 
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